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Abstract. We discuss the processes underlying the excitation of fishbone-like internal 
kink instabilities driven by supra-thermal electrons generated experimentally by 
different means: Electron Cyclotron Resonance Heating (ECRH) and by Lower Hybrid 
(LH) power injection. The peculiarity and interest of exciting these electron fishbones 
by ECRH only or by LH only is also analyzed. Not only the mode stability is explained, 
but also the transition between steady state nonlinear oscillations to bursting (almost 
regular) pulsations, as observed in FTU, is interpreted in terms of the LH power 
input. These results are directly relevant to the investigation of trapped alpha particle 
interactions with low-frequency MHD modes in burning plasmas: in fact, alpha 
particles in reactor relevant conditions are characterized by small dimensionless orbits, 
similarly to electrons; the trapped particle bounce averaged dynamics, meanwhile, 
depends on energy and not mass. 
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1. Introduction and Background 

Fishbone - like internal kink instabilities driven by electrons have been observed for 
the first time on DIII-D in conjunction with Electron Cyclotron Resonance Heating 
(ECRH) on the high field side [1]. There, the excitation was attributed to barely trapped 
supra-thermal electrons, which are characterized by drift-reversal and can destabilize a 
mode propagating in the ion diamagnetic direction in the presence of an inverted spatial 
gradient of the supra-thermal tail. Similar but higher frequency modes were observed 
in Compass-D [2] during ECRH and Lower Hybrid (LH) power injection, with chirping 
frequency comparable with that of the Toroidal Alfven Eigenmode p] (TAE), u ^ ujtae- 
Observations of electron fishbones with ECRH only [U [5] and LH only [6l [7] have been 
also reported in HL-IM and FTU, respectively. More recently, electron fishbones have 
been observed in Tore Supra [8] due to resonant excitation of a double-kink mode by 
supra-thermal electrons generated with LH power injection. 

In the present work, we analyze the peculiar features of electron fishbones versus 
those of the well known ion fishbone [9l [TOl [11]. Due to the frequency gap in the 
low-frequency shear Alfven continuum for modes propagating in the ion diamagnetic 
direction [TT], effective electron fishbone excitation favors conditions characterized by 
supra-thermal electron drift reversal, consistently with experimental observations. For 
the same reason, the spatial gradient inversion of the supra-thermal electron tail is 
necessary, explaining why ECRH excitation is observed with high field side deposition 
only [U m El [121 [13] . Here, we also discuss the pecuhar roles of circulating supra-thermal 
electrons for electron fishbone excitations with LH only: the barely circulating population 
providing directly the mode drive and the well circulating particles controlling both the 
drift-reversal condition as well as the ideal MHD stability via their effect on the plasma 
current profile. The role of LH current drive in controlling sawtooth oscillations via 
the local magnetic shear at the q = 1 surface [q being the safety factor) was recently 
documented by the HT-7 tokamak p3]. 

As in the case of ion fishbones, two branches of the electron fishbone are shown to 
exist: a discrete gap mode [11] and a continuum resonant mode [10]. Contrary to the gap 
mode, the continuum resonant mode can propagate in the electron diamagnetic direction 
as well. Thus, it does not require neither drift-reversal nor inverted spatial gradient of 
the supra-thermal electron tail. However, its threshold condition in this case is higher 
and it requires high power densities to be excited. So, even the case of the continuum 
resonant fishbone mode tends to favor the branch propagating in the ion diamagnetic 
direction, which minimizes continuum damping. If the effective temperature of the 
supra-thermal electron tail is sufficiently high, the present theory predicts that fishbone 
oscillations can be excited at frequencies comparable with those typical of the Geodesic 
Acoustic Mode (GAM) [I5] or the Beta induced Alfven Eigenmode (BAE) [Ml [IT]. 
Unlike the case of fishbone gap modes in the ion diamagnetic gap [11] of the low- 
frequency shear Alfven continuum, fishbone gap modes in the BAE gap [18j do not 
favor the propagation in the ion diamagnetic direction, since the gap structure is nearly 
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Figure 1. Time evolution (from top to bottom) of thermal electron temperature, 
plasma line density, LH coupled power, fast electron temperature fluctuations and 
central radiation temperature in FTU shot # 20865. It is clear that the nonlinear 
behavior of fast electron temperature fluctuations (electron fishbone) reflects the level 
of LH power input. 



symmetric in frequency |19J. Here, we discuss these issues using one single general 
fishbone-like dispersion relation [20l |2T] , describing mode excitation by trapped as well 
as circulating supra-thermal electrons in both monotonic and reversed magnetic shear 
equilibria [22] . 

In this work, we also analyze the nonlinear physics of electron fishbones, of which 
FTU experimental results provide a nice and clear example (see Figure 1): during 
high power LH injection, an evident transition in the electron fishbone signature takes 
place from almost steady state nonlinear oscillations (fixed point) to regular bursting 
behavior (limit cycle). Here, we present a simple yet relevant nonlinear dynamic model 
for predicting and interpreting these observations. 

These results are directly relevant to the investigation of trapped alpha particle 
interactions with low-frequency MHD modes in burning plasmas: in fact, alpha particles 
in reactor relevant conditions are characterized by small dimensionless orbits, similarly 
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to electrons; the trapped particle bounce averaged dynamics, meanwhile, depends 
on energy and not mass. Rigorously speaking, the same argument applies to barely 
circulating particles as well, whose definition is given in Section O For these reasons, 
we could draw a symmetry between trapped ion (alpha particle) and trapped electron 
dynamics and exploit the combined experimental use of LH and ECRH analogously to 
what is done with Neutral Beam Injection and Ion Cyclotron Resonance Heating. 

2. Mode dispersion relations 

The fishbone dispersion relation can be obtained by the standard matching procedure 
of mode structures in the ideal region and inertial layer [23] and generalizing the results 
therein. Here, we choose to solve quasi-neutrality and vorticity equations following the 
procedure of Ref. [24], where the solution of the kinetic layer equations in the Fourier 
space are matched to the ideal region. Letting x = —keir — r^), with ke the poloidal 
wave vector (here, ker^ = —1) and the radius of the q{rs) = 1 surface, we introduce 
the representation 

6(f){x) = j dr]e-'''^6<!>{r]) (1) 

for the scalar potential fiuctuation and other fields. For finite shear, s = rsq'g/qs, the 
asymptotic ideal region solution for odd parity modes when approaching the inertial 
layer is [21] 

\t]\5^ ~ -i{A5<P/27r)iko/\ko\) iri/\r]\) (l + \r]\SW/s^') , (2) 

for 1 <^ \r]\ <ti \'~f/u!\^^, with 7/0; the normalized mode growth rate 
and 6W the normalized potential energy, which, for low-pressure and high 
aspect ratio tokamak equilibria with circular fiux surfaces, is given by SW = 
{2\kg\Ro/B^){6W/\A6^r\^){q^{rs)/rs) [10]. Here, Rq is the tokamak major radius, Bq 
the on axis magnetic field, A6^r is the jump of the radial mode displacement across 
the inertial layer (with A6(j) the corresponding jump in 6(f)) and 6W the mode potential 
energy. Meanwhile, the inertial (kinetic) layer solution is [20] 

\r^\6^ ^ -i{A6<P/27T){kg/\kg\) (t^M) (1 + i\r]\A/\s\) . (3) 

Here, A is the generalized inertia term introduced in [20j. Given Eqs. ([2]) and (j3]), the 
mode dispersion relation reads [TOl [H] 

iA|s| =6W = 6Wf + 6Wk , (4) 

where the fiuid 6Wf, in its simplest expression, is given by ^25] 

6Wf = SnAqo (13/144 - /3%) [r]/ Rl) (5) 

with (3ps = -{Ro/rlf J^' r'^{dp/dr)dr, Ago = 1 - g(r = 0) and /3 = SttP/B^ the ratio of 
kinetic and magnetic pressures. The fiuid term, 6Wf, includes the contribution of the 
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energetic (hot) particle adiabatic and convective responses as well (TU]. Meanwhile, the 
kinetic SWk is [lO] 

5Wk = A^mwl'^r-dr!EdEd\ V i^^MW^ i^W^^^^^ , (6) 
^0 n Jo q J - w 

where m is the energetic particle mass, Uc = {cB/mc) is the cyclotron frequency, 
S = t;2/2, A = fiBo/S = {Bo/B)vl/v^, B ■ V(C - q{r)e) = 0, C is the "toroidal 
angle" chosen such that (r, 6, Q is a toroidal flux coordinate system with straight 
fleld lines (g = g(r)), 1^ = v^^di)-^ j v^\...)di denotes bounce-averagmg, t 
is the arc length along the equilibrium B-fleld, is the bounce/transit time for 
magnetically trapped/circulating particles, Ud is the magnetic drift frequency and 
QFo = {code + u;^)Fo, Cj^Fq = uj-\k x B/B) -VFq, with Fq = Fo{S, \,vii/\v\\) the fast 
particle equilibrium distribution function. In deriving Eq. ([6]), we have closely followed 
[To] and solved for the energetic particle distribution function 

m at m at m uj 

neglecting flnite orbit widths and separating both adiabatic (oc dsFo) as well as 
convective (oc QFq) responses. In this way, one obtains 

SK = ±^£?S!^^„(.)eH<-.W») , (8) 

m UJ uod — LO 

with 5(j) = 6(f)o{f) exp i(C — 6). In this form, the dispersion relation neglects the thermal 
ion kinetic response in the ideal region p6], whose analysis is outside the scope of 
this work and, for our purposes, we can consider as included in the expression of 6Wf 
along with the contribution of the energetic (hot) particle adiabatic and convective 
responses [10]. Note that Eq. ([6]) depends only on the fast particle energy: the only 
residual mass dependence would be through their flnite orbit width, which are neglected 
in the present treatment since we are assuming that the characteristic orbit size is much 
smaller than the fluctuation wavelength in the ideal region. This fact conflrms our 
conjecture that experimental studies of electron flshbones are relevant for understanding 
alpha particle dynamics in burning plasmas, as noted in Section [H More detailed 
discussions on this issue are presented in Section HI 

For the s = case but flnite 5^ = Eq. (j4]) cannot be applied; meanwhile, 

the asymptotic expressions corresponding to Eqs. ([2]) and ([3]) are changed to 

^ -i(A50/2vr)(V|fc.|) iv/\v\) (l " {2kyj3q',M'6W / S') , (9) 
for the ideal region, while the inertial (kinetic) layer solution becomes [27] 

|r/|5$ ~ -i{A5(l)/27T){ke/\ke\) iv/\v\) (l + i«i«2(«i + «2)|r/|V6) , (10) 
where, ai and a2 are deflned as 

«i= -(-^(A + %g.i?o))''' , 

«2= (^f^(A-M^^o))'^' • (11) 
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By asymptotic matching procedure between Eqs. and (fTU]) . we readily derive the 
mode dispersion relation with a simple inertial layer at [22] 



S {Aql - A^)'/' 1 + Ag,/v/Ag2_A2 ^" = sWf + 6Wk 



1/2 



(12) 



with k\\sqsRo = Ag^ = Q's — 1 in this case. 

Equations (jl]) and ( |T2i) are the basis for our linear stability studies of electron- 
fishbones. Their general structure in is known; however, we want to emphasize two novel 
aspects: (i) that Eq. describes the resonant excitation of internal kink fluctuations 
by both trapped as well as barely circulating supra-thermal electron tails; (ii) that the 
analysis of the generalized inertia term, A [20l [21] , demonstrates the existence of ion- and 
electron- fishbones at frequencies comparable with that of GAM [T5] and B AE flGl [T7] . 
That well circulating supra-thermal electron tails can control the internal kink stability 
via their influence on the g-profile, i.e. 6Wf, has been noted for explaining recent 
observations on the HT-7 tokamak [14J and will be simply assumed in this work. 



2.1. Resonant excitation by trapped and barely circulating supra-thermal tails 

For analyzing the different roles of trapped and circulating particles, we move from 
{S, A) to {S, K^) space, with 

1 - (1 - r/i?o)A 

< 1 [0 < A < (1 — r/Ro)] indicating circulating particles, while trapped particles 
have > 1 [(1 — r / Rq) < A < (1 + r / Rq)]. Using the (s, a) model tokamak 
equilibrium [28] (a = —Roq'^dp/dr), the following expressions for the (transit, bounce) 
time of (circulating, trapped) particles are obtained: 

{2Sf^ r 2(r/i?o) 



K 



2{r/Ro) + {l-r/Ro)K^ 



(14) 



AJKiK) ' 41K(l//t) J qRo 
Here, IK stands for the complete elliptic integral of the first kind. In the same way, the 
bounce averaged precession frequency Ud can be computed as [29l [30] : 



S 



{q/r){K^ + 4:r/Ro) 



u,Ro2{r/Ro) + {l-r/Ro)K^ [ 



lE(fi:) 



Ar/R^y 2q 



a 4 



TT 



lE(/t) 



JK(k) 2]K(k] 



(15) 



for circulating particles (/t^ < 1), whereas, for magnetically trapped particles 



8 q 



21E(l/«:) fMlM + i 



U(l/'^) 



K(l/«:; 

^-l/.^ + (2/.^-l)M 



1 



(16) 



UcRor 
a 

3 ' ^'lK(l//t) 

where IE stands for the complete elliptic integral of the second kind. By direct 
inspection of Eqs. ([T^ to ([T^ and accounting for the fact that / d£dX = 
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J d£dK\2r/Ro)[2{r/Ro) + {l-r/Ro)K^]~' by definition of K?, we see that only 
circulating particles with (r/RoY^"^ ^ < 1 contribute to 5Wk on the same footing as 
trapped particles with > 1. Meanwhile, is the strength of the poloidal modulation 
of the parallel velocity along the particle trajectory; thus, we denominate circulating 
particles with (r/RoY^'^ ^ < 1 as barely circulating to distinguish them from the 
well circulating particles with k"^ < {r/R^y/'^. The peculiar roles of trapped and barely 
circulating particles will be further discussed in Section 12.21 in connection with the 
generalized inertia term. A, appearing in Eqs. (jl]) and ( |T2l) . 

Equations (jl]) and (fT2|) generalize the electron fishbone dispersion relations, 
analyzed recently [12], ]13] [31], to both trapped and barely circulating fast particles, 
including (s, a) model equilibrium effects on uj^- A detailed discussion of the circulating 
electron effect on 5Wk was recently given in [32] . A further extension of Eqs. (|4]) and ( fT2l) 
to a broader frequency range than that usually assumed near the ion diamagnetic 
gap |lT] in the low-frequency shear Alfven continuum is discussed in Section 12. 2[ 



2.2. Generalized inertia and high-frequency fishbones 

For the present scope, we need an explicit expression of the generalized inertia term, 
A, appearing in Eqs. (jl]) and (fT2l) . for two limiting cases: (i) the banana regime, 
\uj\ <^ ojhi <^ ujti, with ujbiiojti) the thermal ion bounce (transit) frequency, where [26l 133] 

A^ = {uj^/u\) (1 - uj^pju) [I + (l.6(i?o/r)^/2 + 0.5) q^] ; (17) 

(ii) the high frequency regime, |ci;| ^ uju, where [19] 

ul^E (46/49) + (32/49) (T,/T,) + (8/49) (T^/T,) 2" 



s2_^ ^BAE 

,2 , ,2 



^A 



1 + 



(f^^ (l + (4/7)(Te/T,))' 

Here, oja = VA/{qRo), va is the Alfven speed, u^pi = k x B/i? • VPi/{nimiUJci), k 
is the wave-vector, ujbae = (l^ti{7 / A + Tf./Ti)^^'^ and uju = {2Ti/miY^'^ /{qRo). The 
shear Alfven frequency gap is given by the condition IReA^ < [20] [21] (A is generally 
complex), while the shear Alfven continuous spectrum is described by [IH] 

A^ = qq^Rl . (19) 

The correct form of the enhancement factor oc q^ in Eq. (ITTj) was first pointed out 
in [26]: the l.Q{RQ/rY/'^q^ factor comes from trapped > 1, and barely circulating 
particles, 1 > ^ [r/RoY^"^', the 0.5^^ term, meanwhile, is due to well circulating 
particles, ^ {t/RqY^'^ (see [Appendix A| for a more detailed discussion). It differs 



from the well known 2q^ factor [31] due to the intrinsic limitation of the ideal MHD 
model in assuming an isotropic pressure response: 2q^ would be the result for 5P = 6P\\, 
while 6P± 7^ 6P\\ for the geodesic curvature dynamics in toroidal systems. The problem 
of the kinetic bulk ion inertia enhancement for low frequency (banana-regime) MHD 
modes was analyzed in Refs. [33] [35] [36] , where estimates were given for both inertia 
enhancement as well as ion Landau damping. A more systematic analytic approach 
was given in Refs. [26] and [37]. More recently, it was pointed out that ion Landau 
damping due to the precession resonance with thermal (bulk) ions may be of crucial 
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importance in determining the internal kink mode stability in ITER Here, it is 

worthwhile noticing that the inertia enhancement factor is identical to the zonal flow 
(ZF) polarizability induced by Ion Temperature Gradient (ITG) turbulence [39l HO] . 
This is not a coincidence and is due to the fact that, at long wavelengths, shear 
Alfven wave compressibility due to geodesic curvature coupling at k\\ = is identical 
to the corresponding dynamics of electrostatic waves with = kg = 0, provided that 
diamagnetic effects are neglected. For this reason, we must expect a correspondence 
between ZF polarizability and shear Alfven wave inertia enhancement in the banana 
regime, as in Eq. (fT7|) : a similar correspondence is expected between GAM and Eq. ( !T8|) . 
as pointed out in [21], HI] (see also the following discussion) . 

Similar considerations apply for Eq. f|T8l) . for oju -C |ti;| -C cu^i, where the oc 
term is different from (2g^)~^, predicted by ideal MHD jlSj. It was proposed in |42j . 
within the limits of a local approximation (i.e. without the proof of the existence 
of unstable eigenmodes), that compressibility effects, associated with wave-particle 
resonances due to the periodic toroidal transit motion of thermal ions, may be a source 
of instability for short wavelength shear Alfven waves. Later, other authors |13| HI] 
numerically demonstrated the existence, well below the ideal stability threshold, of 
electromagnetic instabilities due to ion magnetic drift resonances, assuming the very 
short wavelength limit \ujti\ <^ \uj\ \ujdi\- The effect of ion transit resonances was 
reconsidered in ^51 H6j . where it was demonstrated numerically that the uj = uju 
resonance has analogous effects to those of = Udi, and in [17], where these effects 
on resistive interchange modes were analyzed. All these analyses of short wavelength 
drift-type modes are important for the present investigation since the inertial (kinetic) 
layer physics is the same at high and low mode numbers [24j and, therefore, they 
are relevant for the computation of the renormalized plasma inertia for low frequency 
MHD fluctuations. For the same reason, these studies were readily extended to the 
investigation of long- wavelength (low-mode- number) MHD modes [371 HHl HH]- With 
the same expression of A, derived in Refs. [121115], Ref. [19] demonstrated the existence 
conditions of fluctuations of the shear Alfven branch, excited by both energetic as well 
as thermal ions below the ideal MHD stability threshold, based on the general fishbone- 
like dispersion relation [20l [2T] in the form of Eq. (jlj). In the long wavelength limit, 
the expression of A of Refs. [IHl 1121 HS] accounts for the inertia enhancement as well 
as ion Landau damping for uju -C |ci;| <^ uja- For uJu <^ \uj\, it reduces to Eq. (ITSll 
(see also [37]), with an exponentially small ion Landau damping, oc exp(— tu^/tu^^j). 
This favors the formation of fishbone gap modes near the BAE accumulation point for 
conditions with ojbae ^ ^u, i-e. Tg/Tj ^ 1 and/or g ^ 1 [27J. Note that, due to the 
symmetry of the frequency gap described by Eq. (ITS!) . fishbone gap modes near the BAE 
accumulation point can be equally excited in both ion as well as electron diamagnetic 
directions. Meanwhile, the existence condition for the "BAE" -fishbone gap mode is given 
by He {sWf -\- 8Wk^ < [2ni[2I]- That the shear Alfven continuum accumulation point 
(A^ = 0) given by Eq. ( |T8l) is degenerate with the GAM frequency [15], as pointed out 
in [211 IH] , can be verified by direct comparison with the kinetic expression of the GAM 
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frequency given by Ref. [SU]. The degeneracy of BAE accumulation point and GAM 
frequency has been recently noted also in Ref. |51j . 

3. Linear excitation of electron fishbones 

In this Section, we examine more closely the excitation of electron-fishbones on the 
basis of the mode dispersion relations, Eqs. (jl]) and f|T2l) . introduced and analyzed 
in Section [2l We also discuss some experimental evidence of both low- as well as high- 
frequency fishbones, for which the generalized inertia term is given by Eqs. (1171) and ( |T8l) . 
respectively. This frequency classification strictly applies to discrete gap modes, which 
tend to be excited nearby the shear Alfven continuum accumulation points. It can 
be extended to continuum resonant modes as well, when the mode drive is sufficiently 
weak that proximity to accumulation points matters for minimizing continuum damping. 
Generally, strongly driven continuum resonant modes can be excited regardless the shear 
Alfven continuum structure. 

3.1. Low-frequency fishbones 

The crucial features of low-frequency electron fishbone excitations are dictated by the 
asymmetry of the shear Alfven continuum structure at low frequency quantitatively 
expressed by Eq. ( TTTl) . which favors the excitation of modes propagating in the ion 
diamagnetic direction. Consistently with experimental observations [Illll[5], high field 
side ECRH fulfills this requirement and guarantees both drift-reversal of the barely 
trapped supra-thermal electrons as well as the inverted spatial gradient of the supra- 
thermal tail {uj^,/uj > 0) necessary for effective mode excitation. The case of mode 
excitation by LH only |6l [7] follows the same physics with few additional twists. The 
fast electron population which effectively excite the mode are the trapped and barely 
circulating particles (k^ ^ {r/RoY^"^), because of Eqs. (Ej) and (|T5l) . Meanwhile, LH 
power forms a parallel as well as a perpendicular fast electron tail (via Coulomb 
collisions), which is moderately slanted toward the counter-current direction; i.e., despite 
that it guarantees the inverted spatial gradient of the supra-thermal tail {oj^/uo > 0), it 
is less selective than high field side ECRH in producing particles with drift-reversal. In 
the case of mode excitation by LH only [HI [7] , the presence of circulating supra-thermal 
particles is crucial for two reasons: (i) barely circulating particles (k^ ^ {r/RoY^"^) 
effectively contribute to the mode excitation as described by Eq. ([6]); (ii) well circulating 
particles (k^ ^ (r/i?o)^^^) modify the current profile, eventually reversing the magnetic 
shear and broadening the fraction of trapped particles characterized by drift reversal, 
as shown in Eq. f|T6|) . Note that this effect modifies directly the kinetic contribution to 
the internal kink potential energy and is not associated with the MHD (fiuid) potential 
energy change, controlled by LH power via current profile modification, as recently 
discussed for explaining HT-7 observations [H]. As in the case of ion fishbones, two 
branches of the electron fishbone exist: a discrete gap mode [11] and a continuum 
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Figure 2. Absorbed LH power density (broken line) and q profile (solid line) at 
t = 220 ms, as predicted from transport simulations of FTU shot # 20865. The total 
absorbed LH power is Plh = 0.76 MW. 



resonant mode [1^. The latter does not generally require neither drift-reversal nor 
inverted spatial gradient of the supra-thermal tail; however, it has a higher excitation 
threshold and, thus, it is unfavored, particularly for the branch propagating in the a;*e 
direction. 

Applying Eq. ( fT2l) to FTU shot # 20865 (see Figure 1), the almost steady oscillation 
of the mode in the low LH power phase and the absence of sawtooth oscillations suggest 
that 1 ^ Ag^ > 0. This is consistent with the g-profile reconstruction by transport 
simulations, reported in Figure 2 (FTU has no q profile measurements near the magnetic 
axis). Even in the high LH power phase (Figure 3), the minimum-g value remains 
extremely near unity. From experimental observations, ~ 60 krad/s, uj^^i ~ 23 
krad/s, cUbj ~ 70 krad/s, ujh ~ 400 krad/s, ujbae — 900 krad/s and ua — 9.5 Mrad/s. 
Thus u^pi < u ^ Ubi ^ ujti and we can apply Eq. (|T71) . showing > 0. Given the 
oj < ujbi condition, a further generalization of Eq. f|T7|) would be necessary for a rigorous 
analysis including mode damping by precession [38j and precession-bounce resonances 
with thermal ions. These results, however, would simply lead to a redefinition of the 
mode excitation threshold (see Eq. (125|1 below) at the expense of technical complications; 
thus, they will be reported elsewhere. Given Eq. ( fT2l) . for A^ > Ag^ the mode can be 
considered as continuum resonant mode p!Q], following the standard classification |10j . 
Meanwhile, A^ < Ag^ would correspond to a gap mode [H]. In either case we can assume 
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Figure 3. Absorbed LH power density (broken line) and q profile (solid line) at 
t = 280 ms, as predicted from transport simulations of FTU shot # 20865. The total 
absorbed LH power is Plh — 1-69 MW. 



~ Ag^ = 0(10"*^), consistent with g-profile reconstruction by transport simulations 
and with experimental observations. At larger values of Ag^, the mode frequency near 
the accumulation point would rapidly increase up to the BAE frequency and Eq. f|T8l) 
would apply, rather than Eq. f|T71) . Besides the obvious consequence of increasing the 
MHD stability of the system, i.e. 5Wf^ this fact would imply that higher effective 
supra-thermal electron temperature are needed for both balancing 6W j by (51Reiyfc 
and for efficiently driving the mode via wave particle resonances (see Section 13. 2p . 
Altogether, we expect that increasing Ag^ increases the stability of the system, as verified 
experimentally on FTU. 

In the case of the gap mode near the accumulation point [H] , the existence condition 
is bWj + 51Reiyfc < and real mode frequency is given by 

f + IRe^W^feV 

(IReA)^ = ^ql - A ^^^^^ ^ , (20) 

while the growth rate is obtained from 

7 = r [g\m^mk - EmA^ , (21) 

where g = - (bW ^ + IRe^W^fc) / (S^Ag.lReA) and V'^ = dM.eA/du - gd]Re6Wk/duj. 
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For the continuum resonant mode [TD], Eq. (I12p can be written as 



iS (A' - Aql) \Aq, - i (A^ - Aql) '^'1 = 5Wf + SWk . 



(22) 



Assuming Aq 0, for simplicity, the mode dispersion relation becomes 



6Wf + JRe6Wk = (5/V2)A3/2 ^ g , 



(23) 



which determines the mode frequency [TO]; meanwhile, the mode growth rate is defined 
by [52] 




where T = ~{Ro/rs){d]Re6Wk/du - 3S/{2V2)A^/^dA/duj)-\ the effective reso- 



Note that the oc Ph,res term in Eq. fl24|) would change sign for the case of mode excitations 
by fast ions. 

Despite the different structures of Eqs. (120|) and (12T|) with respect to Eqs. (!23ll 
and (Ell), their extension to the nonlinear regime follows the same derivation. For 
this reason, we derive the nonlinear amplitude equations describing the fishbone cycle, 
in Section HI limiting specific applications to the simple case of Eqs. ( !23l) and f l24l) . 
Analogous derivations in other more general cases, included in Eqs. (jl]) and f[T^ . follow 
consequently. 

With FTU shot # 20865 data, S = 0.52 and Ph,c ^ 0.34 x lO"'* at t = 
220 ms, while S = 4.4 and l3h,c ^ 4.9 x 10~^ at t = 280 ms. Lower Hybrid 
power deposition computations provide the supra-thermal electron tail distribution 
function [53], which give Jq" {r /rs)drdrPh,res — 0.85 x 10"^ at t = 220 ms (see Figure 4) 
and Jq^ {r /rs)drdrPh,res — 6.0 x 10""^ at t = 280 ms (see Figure 5). This is consistent 
with mode excitations and the transition observed in Figure 1, when the power level 
is stepped from given Plh = 0.76 MW up to P^h = 1-69 MW. The bursting fishbone 
activity is further discussed in Section HI 

3.2. High-frequency fishbones 

At higher frequencies, u ^ uu, Eq. (fT8|) applies instead of Eq. (ITTI) : thus, the asymmetry 
of the shear Alfven continuous spectrum is lost and modes can equally propagate in 
both ion and electron diamagnetic directions. Equation f|T8|) describes the formation of 
the Beta induced Alfven Eigenmodes (BAE) [16] spectral gap: so, electron fishbones 
propagating in the electron diamagnetic direction and normal pressure profiles could be 
excited. More precisely, high power ECRH experiments with on axis resonance would be 
needed, producing sufficiently high effective supra-thermal electron tail temperatures, 
Th, for the fast particle precession frequency to be of the order of the thermal ion 





(24) 
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Figure 4. Contour plot of the supra-thermal electron tail a.t t — 220 ms, as predicted 
from Fokker-Planck computations of FTU shot # 20865. The total absorbed LH power 
is Plh = 0.76 MW. Velocities are normalized to the core electron thermal speed. The 
radial position is r/a — 0.17 and dashed hues indicate the trapped particle region. 




Figure 5. Contour plot of the supra-thermal electron tail at t — 280 ms, as predicted 
from Fokker-Planck computations of FTU shot # 20865. The total absorbed LH power 
is Plh = 1.69 MW. Velocities are normalized to the core electron thermal speed. The 
radial position is r/a — 0.35 and dashed lines indicate the trapped particle region. 

transit frequency. For the above FTU parameters, this would require Ty^ ^ 100 keV, 
to be compared with the usual values ~ 30 keV, as well as Tg ^ Tj for consistency 
(see Section [2.21) . Obviously, at such high energies of the supra-thermal electron tail, 
relativistic effects can be important and should be included in the expression of |13j . 

The existence condition of gap modes in the BAE frequency gap just below the 
continuum accumulation point is given by He {^Wf -|- 5^4^ < [20l El] , as discussed in 
Section [221 Note that these fishbones, possibly excited below the BAE frequency, could 
be equally excited by ICRH induced fast ions but, in that case, they would propagate in 
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UJ \—r\ / OJii 


T / OJfi 




t 


cJa— n/27r 


'-(-^>t.n7 / 27r 

^pz / ' ' 




.24100E+01 


-.07500E-01 


.lOOOOE+00 


3 keV 


68 kHz 


8.5 kHz 


55 kHz 


.22900E+01 


-.06700E-01 


.15000E+00 


3 keV 


65 kHz 


13 kHz 


55 kHz 


.21700E+01 


-.03200E-01 


.20000E+00 


3 keV 


62 kHz 


17 kHz 


55 kHz 


.22900E+01 


-.10700E-01 


.lOOOOE+00 


4 keV 


75 kHz 


9.9 kHz 


59 kHz 


.21900E+01 


-.09700E-01 


.15000E+00 


4 keV 


72 kHz 


15 kHz 


59 kHz 


.20700E+01 


-.06300E-01 


.20000E+00 


4 keV 


68 kHz 


20 kHz 


59 kHz 



Table 1. Theoretical values of the BAE accumulation point, uja=o from A = [TO] , 
as a function of Tj and w^ni ~ a;^pj/(l + 77^). Fixed parameters are Tg = 6 keV, Rq = 3 
m and rji — 2. Values of ion Landau damping, 7, are also reported. 



the ion diamagnetic direction. The observation of high frequency precessional fishbones 
with ICRH in JET [5^155] can be possibly interpreted as evidence of fishbone excitation 
below the BAE frequency, as predicted by theory. One striking evidence that Eqs. (jlj) 
and ( fTSl) describe these physics is Fig. 11 of [55]. In fact, as the diamagnetic fishbone 
get excited and less free energy is available for the excitation of the precessional fishbone 
(modes are less strongly driven), theory predicts that frequency chirping should decrease 
and the mode frequency should get closer to the accumulation point. However, this is 
evidently not the usual accumulation point at oo^^pi, but rather the accumulation point 
described by Eq. (ITSj) . In fact. Fig. 11 of [55] shows the frequency accumulation at 
about 70 kHz. To test this conjecture, we have computed the BAE accumulation point 
in two ways: (a) via the simplified expression ubae = Q^ti{7/4: + Te/Tj)^/^; and (b) 
via numerical solution of A = 0, with A given by Ref. [19], i.e. including both thermal 
ion transit resonances (for the ion Landau damping evaluation) as well as diamagnetic 
effects (finite u!^:pi). For the JET discharge # 54300 {D plasma with ICRH if-minority 
heating), we have taken Tg = 6 keV, i?o = 3 m and r]i = d\nTi/d\nn = 2, obtaining the 
results reported in Table [TJ Values of Landau damping are typically small. Meanwhile, 
comparisons of theoretical frequencies with the experimental value of ~ 70 kHz suggest 
that a realistic estimate for Tj at the q = 1 surface is Tj ~ 4 keV with 15 kHz < u;^,pj < 20 
kHz, in agreement with experimental observations [55]. The crystal spectrometer for 
this case gives Tj = 2.2 keV, which is a lower bound of the ion temperature at the 
q = I surface and approximately 60% of its value, as suggested by normal experience. 
The good agreement of theoretical predictions with experimental observations confirms 
the sound basis of the proposed interpretation of high frequency precessional fishbones 
observed in JET [HU [55] with ICRH as evidence of fishbone excitation below the BAE 
frequency [27]. The scaling of the BAE accumulation point frequency with T^/Ti can be 
used for diagnostics purposes, similar to the approach proposed in Ref. [56j for Alfven 
Cascades. Actually, the results presented here (see Table [1]) and their dependence on 
diamagnetic effects show that a better evaluation of the accumulation point frequency 
can be obtained by solving A = [I9] rather than using u = ojbae |56j, with the 
additional advantage of computing ion Landau damping. In the case of Alven Cascades, 
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of course, the accumulation point at s = should be evaluated using = fc^g^i^Q [T^ . 
as predicted by Eq. f|T9l) . Note that magnetic shear never enters in the accumulation 
point expression, as expected for local oscillations of the shear Alfven continuum and 
explicitly shown by Eqs. (fTT!) and (fT2|) . 



4. Nonlinear amplitude equation 

FTU experimental results (see Figured]) suggest that the level of LH power input controls 
the transition from nearly steady state to bursting electron fishbone oscillations. Here, 
we want to focus on the bursting regime, where we conjecture that the fishbone is 
a continuum resonant mode [10], described by Eqs. (l23ll and (IMll . on the basis of 
preliminary high time resolution Electron Cyclotron Emission (ECE) measurements, 
which indicate that the frequency of temperature oscillations in the bursting phase 
drops by 40 -r- 50% in 2 -h 3 ms. Better resolved data on the mode frequency chirping 
are needed for more accurate comparisons of theory (this Section) with experiments. 
However, preliminary analyses support our conjecture that the bursting mode phase 
is associated with the excitation of a continuum resonant mode well above marginal 
stability, i.e. |/3h,resmax - Ph,resmm\ ~ f3h,c in Eq. (JJl). For this reason, we expect 
that particle nonlinearities are dominant in dictating the time evolution of the fishbone 
cycle, as recently shown in Ref. [57]. The role of mode- mode couplings on the fishbone 
dynamics was specifically discussed in Ref. [58j . 

For strongly driven fishbone modes, fast particle dynamics is secular in the radial 
direction due to the mode- particle pumping process, originally proposed in Ref. [59] . 
Since, in this case, there is no time for the particles to experience trapping in the 
potential well of the wave, we use a different approach with respect to that of Ref. [60j . 
which postulates proximity to marginal stability and describes the nonlinear evolutions 
of modes with slowly varying frequencies due to structures in phase space near particle 
resonances. Here, we adopt the 4-wave modulation interaction model, introduced 
by Chen et al [61] for analyzing modulational instabilities of the radial envelope 
of Ion Temperature Gradient driven modes in toroidal geometry, extending it to the 
modulations on the fast particle distribution function due to nonlinear mode dynamics, 
as proposed in Ref. [62]. In the following, we show that the resonant particle motion 
is secular with a time-scale inversely proportional to the mode amplitude. In order 
to qualitatively compare the model predictions with FTU experimental results on the 
fishbone repetition rate, we show that our nonlinear model equations are expressible 
in terms of a predator-prey like model with a limit cycle. This model differs from the 
existing qualitative models (Refs. [10] and [TT]) in that it is structurally stable, i.e., the 
periodic dynamics not destroyed in the presence of higher-order perturbation terms (see 



Appendix B). The transition to the stable limit cycle behavior occurs via a marginal 
oscillatory regime (i.e., the center ii to use the proper terminology), which is structurally 
unstable and is also revealed in those models discussed in Refs. jlO] and [TT]. 



We can generalize Eqs. (jl]) and (fT2i) to include supra-thermal electron tail nonlinear 
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dynamics by closely following the procedure of Ref. [53] • In the present treatment, as 
discussed above, we choose to neglect fishbone nonlinear dynamics associated with mode- 
mode couplings. For the case of continuum resonant fishbones [lO], this approximation 
allows us to retain the fundamental dynamics [211 EZ] to make significant analytic 
progress, as shown below. 

Under the action of the fishbone mode, the toroidally and poloidally symmetric 
(zonal [SHUT]) nonlinear modification of the fast electron distribution function, Eq. ([7]), 
can be obtained from the nonlinear gyrokinetic equation [M] and is given by [63] 



V i—k' — 



d 



dr 









6(f)k'6Hk" 







(26) 



k^=k'+k" 

where k'^ = —k'^, kg = —kg, we have neglected finite electron orbit widths and assumed 
SE\\k = 0. By direct substitution and using Eqs. ([7]) and ([H]), Eq. (1251) is readily reduced 
to 

d 
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SH, 



NL,: 



2 2 d 
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.(27) 



Here, v y exp i(l — g)6' = for trapped particles and = ^iro/fs is the normalized 
radial displacement of the mode, which is assumed to be the usual step function. The 
presence of the imaginary part of the particle response on the RHS of Eq. ( |271) indicates 
the crucial roles played by resonant particles [65]. Meanwhile, by definition of the QFq 
operator {QFq/u ~ d^Fo + kg/ {ujuj^drFo), the RHS contains both oc drF^ and O^Fq 
terms [65]. Thus, integrating both sides in velocity space, Eq. (127|) can be easily put 
in the form of a diffusion equation describing the relaxation of the fast particle profile 
within the g = 1 surface: 



d_ 
dt 



nh = Nh ujcuj" 
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r r„ 
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Here, Nh indicates the fast electron source term due to additional power input, we have 
defined the effective fraction of fast electrons feff,h and 
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Qres^h 



UJ 



(Fo) 



-1 



k\\V\\ 



UJ 



Em 



'QFo' 



UJ 



(29) 



having indicated velocity space integration by angular brackets. From Eqs. fl28l) and fl29l) 
we recognize that the nonlinear diffusion coefficient due to the fishbone within the q = 1 
surface is given by Dnl ^ 2ur'^Jeff,h \S^of. 

One obvious consequence of Eq. (127|) is the time evolution of the supra-thermal 
electron tail contribution to 6Wf via their convective responses. In fact, nonlinearly [63] 



QFq QFq H --;t-6HnL,z 

uJr or 



(30) 



Thus, the expression of dtSWf^^L is readily obtained from that of 6Wf by direct 
substitution of the supra-thermal electron tail distribution, drFo, with the expression of 
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dtdrSH]^L z from Eq. (1271) . For this reason, in the present work we will simply assume 
it as given, without providing further detailed discussions. The other effect of Eq. ([2] 
is to introduce a nonlinear modification to Eq. ([8]) in the form 



6K 



NL 



i{C-g(r)e) 



(31) 



Bq uod — oj ijJ dr 
Using Eq. (!27|) . meanwhile, the nonlinear modification for the resonant contribution 
(imaginary part) of 5Wk is obtained in the form: 



d_ 
di 
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k,NL 
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(32) 



Note that here, as in the following Eqs. (!35|) . (!37|) and (1391) . partial time derivation is 
intended at constant frequency, which, in the present problem, can be itself a function 
of time and vary on the nonlinear time scale. An intuitive derivation of Eq. fl32|) can be 
obtained from Eqs. ([6]), ([8]) and fl3T]) . noting that, for resonant particles involved in the 
6Hnl,z dynamics 

i 



(ud - uj) SK 
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Recalling the definition of drPh,res, given below Eq. 1^^, i.e. 
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Eq. ( l32l) can be approximated as 
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(34) 



(35) 



Here, for simplicity, we have assumed that the radial variation of uJd is essentially oc (1/t) 
within the q = 1 (or minimum-g) surface; meanwhile, C is a constant which may be 
computed exactly, given Eq. (!32|) and the definition of Ph,res, Eq- (|34l) . Using these 
results, the amplitude evolution equation can be formally written as Eq. ([2i 



(d/dt) 15^0 1^ = 2r 



(r/r J idph,res/dr) dr - ph,, 



(36) 



where the nonlinear evolution equation for the resonant fast particle pressure gradient 
becomes 
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Here, Ug^t is the reconstruction rate of Ph,res, i-e. 

vr^ mujc r'^k, 




''ii/l'"iil=±i 

where dtFo^^xt is the rate of change of the fast particle distribution function due to 
external sources (inclusive of Coulomb collisions). 

Following the same formal steps adopted for the derivation of Eq. ( [32l) . we can 
obtain the expression of M,eSWk,NL at the next order in the asymptotic expansion in 
{■j/ujI ^ l/\u!TiyL\, whcrc Tjyi is the nonlinear time scale: 

r2 D J2. 



^lRe(5Wfc ~ — 2^mu;fu;^ l —dr [ dSdX \^ —e 
9t ' ri Jo q J ,f-'_^^ uJd 



or or 



D||/|l>|||=±l 



^1/2 d£ 



X ei^Me^^^e-'^QFo Im rV |5^o| 



1 



(39) 



The real frequency of the fishbone mode in the nonlinear regime is still given by an 
equation in the form of Eq. fl23l) : i.e., the mode frequency is expected to chirp downward 
as the fast particles relax, according to Eq. (1251) . More specifically, the nonlinear 
evolution equation for the real frequency is 

bWf + (5#/,7VL + lRe(5W^fc + IRe^W^fc,^^ = (5/V2)A^/2 ~ . (40) 

The effect of external sources (inclusive of Coulomb collisions), can be straightforwardly 
added into Eq. fHOj) via time varying lRe(5VFfe and 5W j due to dtFo^exti similarly to 
Eqs. ( !37|) and (!38|) above. Detailed discussions of these issues will be given elsewhere. 

That 5Wk,NL is predominantly imaginary, as it emerges from comparisons of 
Eqs. fl32l) and fl39|) . suggests that the nonlinear fishbone cycle is essentially determined 
by the fast particle scattering out of the resonant region. Given the fluctuation 
level of the mode, the nonlinear time scale, t/vl, as derived from Eq. (137|) . scales as 
{2CY/'^ujtnl ~ 1(5^01"^ 5 consistent with the predator-prey model for the fishbone cycle 
proposed in [in] and in contrast with the time scale oc l^^oT^ of Eq. (|39ll . Meanwhile, 
the oc l^^oT^ scaling of characteristic times pTlj is consistent with the time behavior of 
Eqs. ( I28l) and (HOl) . describing, respectively, the fast particle relaxation and the time- 
dependent nonlinear frequency shift of the mode. This picture of the fishbone nonlinear 
dynamics agrees well with the mode-particle pumping process, originally proposed in 
Ref. [59] . 

Equations fl36l) to fHOl) describe the bursting fishbone cycle when the dynamics is 
dominated by coherent nonlinear interactions, typical of the mode-particle pumping 
process [59], in the presence of quasilinear wave-particle resonances. This formal 
analysis, thus, is equivalent in spirit to the approach of [66] and the numerical analysis 
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of [HZ], but has the advantage of treating exphcitly the energetic particle nonhnear 
dynamics. They can be analyzed with different levels of complexity and their detailed 
analyses will be reported elsewhere, along with comparisons with FTU experimental 
observations. Here, we want to emphasize that Eqs. (l36l) and (!37|) are already a 
simple yet relevant model which describes the fishbone cycle when the dynamics due 



to nonlinear frequency shifts is neglected [TOl [TT]. In [Appendix B[ we show that 
these equations can be reduced to a predator-prey system, similar to the ad-hoc model 
introduced in Ref . [lO] . The main results of that analysis are that the nonlinear system 
is characterized by small oscillation about a fixed point. For increasing LH power 

1 /2 

input, the system approaches a limit cycle of period tff, ~ 27i/ (2TPh,cJ^ext) , with 
Sf3h/Ph,c ~ '^ext /(2r/5?i,c)^''^ estimating the loss of fast particle in one fishbone burst. 
Given these results, the present estimate of tfb is consistent with that of Ref. |10j . 
tfh ^ {SPh/ Ph,c)i^ext- With the parameters corresponding to the high power phase, 
Plh = 1-69 MW, of FTU shot # 20865, we obtain tfh ~ 5-^10 ms. The good agreement 
we obtain on the estimate of the fishbone period (no measurements are available of the 
losses in the perpendicular supra-thermal electron tail) motivates further experimental 
investigations for more detailed comparisons with theoretical model predictions. 

As in Section [31 we may discuss our conjecture of the relevance of electron fishbone 
experimental studies for gaining insights into linear and nonlinear burning plasma 
dynamics. As stated already, the symmetry breaking between fast electron and fast 
ion bounce averaged dynamics is caused by finite orbit width effects (linear dynamics). 
In this respect, the typically small dimensionless orbits of fast electrons may generate 
behaviors analogous to those of well confined fast ions in thermonuclear plasmas. 
When analyzing nonlinear fast particle behaviors, another source of symmetry breaking 
between bounce-averaged fast electron and ion dynamics emerges from Eqs. fl271) . fl32l) 
and fl39l) : the term oc v\\expi{l — q)6 {v\\ exp i{l ~ q)9 = for trapped particles), is 
responsible for the barely circulating fast particle radial transport due to the magnetic 
component of the fluctuations and clearly depends on the particle mass. For particle 
distribution functions that are symmetric in v\\, this term is unimportant. Clearly it 
is not so for the LH driven fast electron distribution, producing a perpendicular fast 
particle tail, which is moderately slanted toward the counter-current direction. This 
physics is embedded in the C factor, introduced in Eq. (l35l) : experimentally varying the 
power mix of LH and ECRH can, thus, control both the excitation condition of electron 
fishbones as well as the amount of radial transport due to the magnetic component of 
the fluctuation. 



5. Discussions and conclusions 

In this work, we have analyzed the excitation of electron fishbones by both trapped as 
well as barely circulating supra-thermal electrons, providing a unified explanation of the 
various experimental observation of these modes. In particular, we have analyzed the 
peculiarities of electron fishbone excitation on FTU by LH power only, explaining the 
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different roles of trapped and circulating supra-thermal electron tails. The possibility 
of exciting fishbone modes at frequencies just below the BAE accumulation point by 
both fast electrons and ions is also discussed and conjectured as interpretation of the 
experimental evidence of ICRH driven fishbone modes in JET, accumulating at finite 
frequency above the diamagnetic gap as the mode drive is weakened. 

We have derived nonlinear amplitude equations, which describe the nonlinear 
bursting fishbone cycle due to the mode-particle pumping process in the presence of 
quasilinear wave-particle resonances when mode-mode couplings are neglected. These 
equations are qualitatively equivalent to a predator-prey like system, whose predictions 
are consistent with the corresponding ad-hoc model, originally proposed for explaining 
the ion fishbone cycle. 

The most interesting feature of electron fishbones is their relevance to burning 
plasmas. In fact, unlike fast ions in present day experiments, fast electrons are 
characterized by small orbits, which do not introduce additional complications in 
the physics due to nonlocal behaviors, similarly to alpha particles in reactor relevant 
conditions. Meanwhile, the bounce averaged dynamics of both trapped as well as barely 
circulating electrons depends on energy (not mass): thus, their effect on low frequency 
MHD modes can be used to simulate/analyze the analogous effect of charged fusion 
products. Symmetry breaking between fast electron and ion bounce averaged dynamics 
is caused by finite orbit width effects (linear dynamics) and by radial transports due to 
the magnetic component of the ffuctuations (nonlinear dynamics of the barely circulating 
particles). In this respect, the combined use of ECRH and LH provide extremely flexible 
tools to investigate various nonlinear behaviors, of which FTU experimental results 
provide a nice and clear example (see Figure 1). 
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Appendix A. Generalized inertia and the peculiar roles of trapped and 
circulating particles 



Here, we further discuss the peculiar roles of trapped and circulating particles in 
determining the generalized plasma inertia, continuing the analyses of Section 12. 2[ For 
this scope, we use the analogy between ZF polarizability and shear Alfven wave inertia 
enhancement in the banana regime, as in Eq. (|T71) . Closely following Refs. [39l HQ], we 
can identify the different roles of trapped and barely circulating particles in determining 
the l.G^Ro/rY^'^q'^ factor in Eq. f|T7|) [26] . It is easily shown that the inertia enhancement 
can be written in compact form as: 
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Here, the lowest order solution SK^'^^ in the Icu/cUf,! asymptotic expansion is |26j : 
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with f II the fluctuating component of the parallel velocity, defined such v\\ 
Eqs. f lA.2p to flA.4p . we readily obtain [26] 
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Here, 6 = 0[(r//2o)^''^]; thus, the first two terms on the right hand side (RHS) represent 
the contribution of barely circulating particles (~ 0.43), while the last term on the RHS 
(~ 1.20) comes from trapped particles. Note that the structure of Eq. ( 1A.6I) is the same 
as that involved in the ZF polarizability [391 HQ], expected. 

At low frequency, |a;| ^ Ubi ~ {r / RoY^'^cou, ion Landau damping is strongly 
decreased due to particle trapping [37j; meanwhile, the inertia enhancement due to 
(well) circulating particles reduces to the 0.5g^ factor in Eq. (fTTj) . Note that, as pointed 
out in Section [2^21 the origin of this factor is crucially related to the anisotropic pressure 
response due to geodesic curvature couplings. In fact, neglecting diamagnetic frequency 
for simplicity, the usual inertia term oc kj{uj^/uj\)6(l){r) is changed into (see Eq. (1A.2I) ) 

kl^6<l>ir) + kl^q' Up\\, + 5P^^ , (A.7) 

where (5P|| j and SP±i are the normalized amplitudes of the oc sin 6 thermal ion parallel 
and perpendicular pressure perturbations due to geodesic curvature. For \u\ <ti uu one 
easily finds [19] 

SP = ^50(r) , 

5Pu = - ^H{r) . (A.8) 



Thus, the 0.5g^ factor is obtained because 6P±i ^ SP±i, while assuming ( 6P\\i + 6P± 
26P\\i = 26Pi would give the usual 2q'^ factor 



Appendix B. A predator-prey like model for the fishbone cycle 

Predator-prey models for the fishbone cycle are known since the original works on the 
resonant continuum [TU] and discrete gap [TT] fishbone modes. Here, we demonstrate 
that the nonlinear model equations for the fishbone cycle, Eqs. ( !36|) to (HOi) . can be 
reduced to a structurally stable predator-prey like model with a stable limit-cycle 
behavior. This gives us a qualitative picture of the fishbone dynamics, which agrees with 
experimental observations on FTU and with earlier work in Ref. [10]. The derivation of 
a structurally stable model for the fishbone cycle offers the mathematical and physical 
foundation of previously proposed phenomenological models [TOl [11] . 

For the sake of simplicity, we neglect the dynamics due to time-dependent frequency 
shifts, contained in Eq. pUj) . and we also assume that Ph,res can be described by a 
characteristic spatial scale A < r^, so that the first radial derivative drPh,res — B/A 
and the second radial derivative dl(3h,res — with E a characteristic value of 

Ph,res- Note that the second radial derivative is set negative in order to guarantee the 
stability of the absorbed LH power density profile. We take for granted that the main 
contribution to the integral in Eq. fl36|) comes from a shell between — A and r^. As 
an illustration, the profile in Figure 3 suggests that r^/a — A/a ~ 0.2 with a negative 
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second derivative onwards. Integrating in Eq. (35) from — A to we find 

^ = -2T{(3-B)A , (B.l) 

where we preferred simpler notations (3 = (3h,c and A = |5^oP- Note that the time 
derivative in Eq. ( ]B.1|) changes sign when B crosses the (3 value. If we now turn to 
Eq. (1371) and approximate the radial derivatives of the profile function Ph,res by their 
characteristic values through the distance A, we get 

(B,2) 

with q = 2Cuj'^TNL{rs/ ■ Differentiating on the left of Eq. ( 1B.2I) we find 

^(^- B + A^(^- uext] B = -^A'B . (B.3) 
dt \dt J dt \dt J tnl 

We are interested in the asymptotic dynamics when the time t +oo and we 
bypass the intermediate transitional-type patterns of behavior when the system basically 
evolves from a linear starting regime into a strongly nonlinear stage when the coupling 
between the dynamical parameters comes into play. Keeping first time-derivatives and 
suppressing higher-order differential terms we rewrite Eq. ( IB. 31) as 

^(^-iye.t]B + ^(^-u,^t]B=-^A''B , (B.4) 

where I/tjvl stands for d/dt where appropriate to constitute the right ordering. In 
writing Eq. ( 1B.4I) . we considered that the dynamical time scale is of the same order 
of the nonlinear time tjvl so that the nonlinearities are essentially present through the 
dynamics. Note that Eq. ( IB. 41) does not admit the linear limit as particular case. We 



also took into account that the time derivative dA/dt is such as to satisfy the dynamic 
Eq. ( IB.lj) . Substituting dA/dt from Eq. ( IB.lj) . after simple algebra one obtains 

'^ = uB-eB^- qAB , (B.5) 

with u = Vext+I3Q and = 2T [1 ~ Vext'TN l) ■ Equations ( 1B.1I) and ( IB. 51) form a predator- 
prey system of equations. If we change the notations in Eqs. (IB. ID and (IB. 51) such that 
^ is X and B is y and introduce the parameters /i = 2(3T and k = 2r, we can represent 
our predator-prey model in the canonical form 

X = —fix + kxy , (B.6) 

y = iyy-Qy'^- qxy , (B.7) 

where the dot denotes time differentiation. Equations (IB. 60 and (\B.7\i define a dynamical 
system with an unstable hyperbolic point at the origin and an equilibrium (i.e., fixed 
point) at xo = u / q — Qfi/qk and yo = fi/k. The unstable hyperbohc point at the origin, 
{x, y) = (0, 0), is a signature of the linear instability of the system. Meanwhile, the fixed 
point at (x, y) = (xq, yo) dictates the nonlinear behaviors. The term with y"^ in Eq. ( 1B.7P 
is important as it guarantees the structural stability of the model, in the topological 
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sense. If O > 0, the system shows a stable hmit cycle behavior, in agreement with the 
bursting mode signatures in FTU, as discussed below. A transition to the limit-cycle 
dynamics corresponds to a nonlinear time tnl, which is comparable to or shorter than 
u'^f-. Physically, this means that the system accommodates the external changes in the 
profile function and, thus, kind of digests the external power density input. If one wishes 
to obtain a marginal estimate, then the procedure is to let be of the order of u'^^ 
and neglect the term with G in Eq. (IB.7p . yielding 

X = —fix + kxy , (B.8) 

y = iyy-qxy , (B.9) 

with the fixed point at Xq = u/q and yo = fi/k. A perturbation analysis of the 
reduced Eqs. ( IB. 81) and (IB. 91) shows that the dynamics are periodic, with frequency 
= (/iz/)-*^/^ = {2TPh^c^exty^'^- The trajectories of the system in the {x,y) phase space 
are defined through |68j 

x'/^y'/'' exp (^-^^±^^ = const (B.IO) 

and are plotted numerically in Figures 6 to 8 for different values of the ratio a = fi/u = 
{2TPh,c/^ext) and with normalized axis scales {x/xo,y/yo): in this way the fixed point 
is always (1,1). Given Sx and Sy the excursions about the fixed point with frequency 
fl, we typically have Sy/y^ ~ z/e^t^^"^ for 5x/xq ~ 1. This means that the characteristic 
excursion of jSh about the fixed-point value is 5Ph/Ph,c ~ {^ext/'2'TPh,cY^'^- With these 
estimates, the period of fishbone burst is tfb ~ 2n/Q ~ 27i{SPh/ Ph,c)^exty in agreement 
with the estimate given in Ref. ^U\. Note that the wider the oscillation amplitude 
around the fixed point is, the more important the non-harmonic behavior becomes in 
the periodic motion of the system, as it is clearly visible in Figures 6 to 8, consistently 
with the electron fishbone burst signature of Figure 1. The nonlinear excursions of the 
system about the fixed point have an amplitude which is dictated by the external power 
density input into the wave-particle resonance region, i.e., ceteris paribus, by Uext, which 
is experimentally controlled via the additional power level. 
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